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ABSTRACT 

We present a calculation of the Sun’s motion through the Milky Way Galaxy 
over the last 500 million years. The integration is based upon estimates of the 
Sun’s current position and speed from measurements with Hipparcos and upon 
a realistic model for the Galactic gravitational potential. We estimate the times 
of the Sun’s past spiral arm crossings for a range in assumed values of the spiral 
pattern angular speed. We hnd that for a difference between the mean solar 
and pattern speed of Gq — Gp = 11.9 ± 0.7 km s“^ kpc“^ the Sun has traversed 
four spiral arms at times that appear to correspond well with long duration cold 
periods on Earth. This supports the idea that extended exposure to the higher 
cosmic ray flux associated with spiral arms can lead to increased cloud cover and 
long ice age epochs on Earth. 

Subject headings: Sun: general — Earth — cosmic rays — Galaxy: kinematics 
and dynamics 


1. Introduction 

Since its birth the Sun has made about 20 cycles around the Galaxy, and during this 
time the Sun has made many passages through the spiral arms of the disk. There is a growing 
interest in determining how these passages may have affected Earth’s environment. Shaviv 
(2002, 2003) makes a persuasive argument that there is a correlation between extended 
cold periods on Earth and Earth’s exposure to a varying cosmic ray flux (GRF). Shaviv 
proposes that the GRF varies as the Sun moves through Galactic spiral arms, regions with 
enhanced star formation and supernova rates that create more intense exposure to cosmic 
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rays. The CRF experienced by Earth may affect the atmospheric ionization rate and, in 
turn, the formation of charged aerosols that promote cloud condensation nuclei (Harrison 
& Aplin 2001; Eichkorn et ah 2002). Marsh & Svensmark (2000) show that there is a close 
correlation between the CRF and low altitude cloud cover over a 15 year time span. Thus, 
we might expect that extended periods of high CRF lead to increased cloud cover and surface 
cooling that result in long term (Myr) ice ages. Spiral arm transits may affect Earth in other 
ways as well. Yeghikyan & Fahr (2004) suggest that during some spiral passages the Earth 
may encounter interstellar clouds of sufficient density to alter the chemistry of the upper 
atmosphere and trigger an ice age of relatively long duration. The higher stellar density 
in the arms may more effectively perturb the Oort cloud of comets and lead to a greater 
chance of large impacts on Earth, and this combined with the possible lethal effects of nearby 
supernova explosions could cause mass extinctions during during passages through the spiral 
arms (Leitch & Vasisht 1998). On the other hand, the record of terrestrial impact craters 
suggests a variation on a time scale shorter than the interarm crossing time, but possibly 
related to the Sun’s oscillations above and below the disk plane (Stothers 1998). 

A comparison of the geological record of temperature variations with estimates of the 
Sun’s position relative to the spiral arms of the Galaxy is difficult for a number of reasons. 
First, our location within the disk makes it hard to discern the spiral structure of the Galaxy, 
particularly in more distant regions. Nevertheless, there is now good evidence that a four-arm 
spiral pattern is successful in explaining the emissions from the star-forming complexes of 
the Galaxy (Russeil 2003). Second, the angular rotation speed of the Galactic spiral pattern 
is still poorly known, with estimates ranging from 11.5 (Gordon 1978) to 30 km s“^ kpc“^ 
(Fernandez, Figueras, & Torra 2001) (see reviews in Shaviv 2003, Bissantz, Englmaier, & 
Gerhard 2003, and Martos et al. 2004). Finally, the Sun’s orbit in the Galaxy is not circular, 
and we need to account for the Sun’s variation in distance from Galactic center and in orbital 
speed to make an accurate estimate of the Sun’s position in the past. 

Here we present such a calculation of the Sun’s path through the Galaxy over the last 
500 Myr. It is based upon the Sun’s current motion relative to the local standard of rest 
as determined from parallaxes and proper motions from the Hipparcos Satellite (Dehnen & 
Binney 1998a) and on a realistic model of the Galactic gravitational potential (Dehnen & 
Binney 1998b). We discuss how the spiral pattern speed is critical to the estimates of the 
times of passage through the spiral arms, and we show a plausible example that is consistent 
with the occurrence of ice ages during spiral arm crossings. 
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2. Integration of the Sun’s Motion 

An integration of the Sun’s motion was made using a cylindrical coordinate system for 
the Galaxy of {R, 0, Z). We hrst determined the position and resolved velocity components 
of the Sun in this system using the velocity of the Sun with respect to the local standard of 
rest (Dehnen & Binney 1998a) and the Sun’s position relative to the plane (Holmberg, Flynn, 
& Lindegren 1997). We then performed integrations backward in time using a fourth-order 
Runge-Kutta method and a model for the Galactic potential from Dehnen & Binney (1998b). 
We adopted the model (#2) from Dehnen & Binney (1998b) that uses a Galactocentric 
distance of Ro = 8.0 kpc and a disk stellar density exponential scale length of = 2.4 kpc. 
This model has a circular velocity at Ro = 8.0 kpc of 217.4 km s“^. We used time steps of 0.01 
Myr over a time span of 500 Myr. Note that the model potential is axisymmetric and does 
not account for the minor variations in the held near spiral arms. We also ignore accelerations 
due to encounters with giant molecular clouds, since their effect is small over periods less 
than 1 Gyr (at least in a statistical sense; Jenkins 1992). The full set of coordinates as a 
function of time is not included here, but interested readers can obtain the digital data from 
our web site^. 

The Sun’s journey in cylindrical coordinates is illustrated in Figure 1. The top panel 
shows the temporal variation in distance from Galactic center, and we see the radial oscil¬ 
lation that is expected from the “epicycle approximation” for nearly circular orbits (Bin¬ 
ney & Tremaine 1987). The period is 170 Myr and the corresponding frequency is 36.9 
km s“^ kpc“^, which is close to the expected value of 36.7 ± 2.4 km s“^ kpc“^ based upon 
the local Oort constants (Feast & Whitelock 1997). The middle panel shows the advance 
in azimuthal position with the orbit (small departures from linearity rehect speed varia¬ 
tions that conserve angular momentum). The Sun has completed just over two circuits of 
the Galaxy over this time span. The lower panel shows the oscillations above and below the 
Galactic plane. The period is approximately 63.6 Myr, but there are cycle to cycle variations 
caused by the varying radial density in the model. This period P is approximately related 
to the mid-plane density at the average radius, p = (26.43 Myr/P)^ Mq pc“^ (Binney & 
Tremaine 1987). The period for our model of the solar motion corresponds to a mid-plane 
density of 0.17 Mq pc“^, which is close to current estimates of the Oort limit of 0.15 ± 0.01 
Mq pc“^ (Stothers 1998). Thus, while the estimates of motion of the Sun in the Z direc¬ 
tion are secure for the recent past, probable errors in the period of approximately 7% may 
accumulate to as much as half a cycle error in the timing of the oscillations 500 Myr ago. 
The errors in the estimates of the Sun’s current Galactic motions (Dehnen & Binney 1998a) 


^http: //www.chara.gsu.edu/"^gies / solarmotion.dat 
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have only a minor impact on these trajectories. For example, the error in the V component 
of motion amounts to a difference of only 3° in 0 over this 500 Myr time span. 

We next consider the motion of the Sun in the plane of the Galaxy relative to the spiral 
arm pattern. The disk of the Galaxy from the solar circle out-wards appears to display a 
four-arm spiral structure as seen in the emission of atomic hydrogen (Blitz, Fich, & Kulkarni 
1983) and molecular GO (Dame, Hartmann, & Thaddeus 2001) and in the distribution of 
star forming regions (Russeil 2003). We show in Figure 2 the appearance of the Galactic 
spiral arm patterns based on the model of Wainscoat et al. (1992) but with some revisions 
introduced by Gordes & Lazio (2003)^. This representation is very similar to the pattern 
adopted by Russeil (2003). We have rescaled the pattern from a solar Galactocentric radius of 
8.5 kpc to a value of 8.0 kpc for consistency with our model of Galactic potential from Dehnen 
& Binney (1998b). Each arm is plotted with an assumed width of 0.75 kpc (Wainscoat et al. 
1992) and each is named in accordance with the scheme of Russeil (2003). The dotted line 
through the center of the Galaxy indicates the current location of the central bar according 
to Bissantz et al. (2003). The pattern speed of the bar may be similar to that of the arms 
(Ibata & Gilmore 1995) or it may be faster than that of the arms (Bissantz et al. 2003), in 
which case the bar - arm relative orientation will be different in the past. 

The placement of the Sun’s trajectory in this diagram depends critically on the relative 
angular pattern speeds of the Sun and the spiral arms. The mean advance in azimuth in our 
model of the Sun’s motion corresponds to a solar angular motion of Dq = 26.3 km s“^ kpc“^. 
If the difference in the solar and spiral arm pattern speeds, Dq —Dp, is greater than zero, then 
the Sun overtakes the spiral pattern and progresses in a clockwise direction in our depiction 
of the Galactic plane. Unfortunately, the spiral pattern speed is not well established and 
may in fact be different in the inner and outer parts of the Galaxy (Shaviv 2003). Several 
recent studies (Amaral & Lepine 1997; Bissantz et al. 2003; Martos et al. 2004) advocate 
a spiral pattern speed of Dp = 20 ± 5 km s“^ kpc“^, and we show in Figure 2 the Sun’s 
trajectory projected onto the plane for this value (Dq — Dp = 6.3 km s“^ kpc“^). Diamonds 
along the Sun’s track indicate its placement at intervals of 100 Myr. We see that for this 
assumed pattern speed the Sun has passed through only two arms over the last 500 Myr. 
However, if we assume a lower but still acceptable pattern speed of Dp = 14.4 km s“^ kpc“^ 
(shown in Fig. 3 for Dq — Dp = 11.9 km s“^ kpc“^), then the Sun has crossed four spiral arms 
in the past 500 Myr and has nearly completed a full rotation ahead of the spiral pattern. 
Thus, the choice of the spiral pattern speed dramatically influences any conclusions about 
the number and timing of Sun’s passages through the spiral arms over this time interval. 


^http://astrosun2.astro.Cornell.edu/'^cordes/NE2001/ 
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The duration of a coherent spiral pattern is an open question, but there is some evidence 
that long-lived spiral patterns may be more prevalent in galaxies with a central bar. For 
example, numerical simulations of the evolution of barred spirals by Rautiainen & Salo (1999) 
suggest that spiral patterns may last several gigayears. Their work suggests that the shortest 
time scale for the appearance or disappearance of a spiral arm is about 1 Gyr. Therefore, it 
is reasonable to assume that the present day spiral structure has probably been more or less 
intact over the last 500 Myr (at least in the region of the solar circle). 


3. Discussion 

Shaviv (2003) argues that the Earth has experienced four large scale cycles in the CRF 
over the last 500 Myr (with similar cycle times back to 1 Gyr before the present). Shaviv 
shows that the GRF exposure ages of iron meteorites indicate a periodicity of 143 ± 10 Myr 
in the GRF rate. Since the cosmic ray production is related to supernovae and since Type II 
supernovae will be more prevalent in the young star forming regions of the spiral arms, 
Shaviv suggests that the periodicity corresponds to the mean time between arm crossings 
(so that Earth has made four arm crossings over the last 500 Myr). Shaviv (2003) and Shaviv 
& Veizer (2003) show how the epochs of enhanced GRF are associated with cold periods on 
Earth. The geological record of climate-sensitive sedimentary layers (glacial deposits) and 
the paleolatitudinal distribution of ice rafted debris (Frakes, Francis, & Syktus 1992; Growell 
1999) indicate that the Earth has experienced periods of extended cold (“icehouses”) and 
hot temperatures (“greenhouses”) lasting tens of million years (Frakes et al. 1992). The long 
periods of cold may be punctuated by much more rapid episodes of ice age advances and 
declines (Imbrie et al. 1992). The climate variations indicated by the geological evidence 
of glaciation are conhrmed by measurements of ancient tropical sea temperatures through 
oxygen isotope levels in biochemical sediments (Veizer et al. 2000). All of these studies lead 
to a generally coherent picture in which four periods of extended cold have occurred over the 
last 500 Myr, and the midpoints of these ice age epochs (lAE) are summarized in Table 1 
(see Shaviv 2003). The icehouse times according to Frakes et al. (1992) are indicated by the 
thick line segments in each of Figures 1, 2, and 3. 

If these lAE do correspond to the Sun’s passages through spiral arms, then it is worth¬ 
while considering what spiral pattern speeds lead to crossing times during ice ages. We 
calculated the crossing times for a grid of assumed values of Gq — Gp and found the value 
that minimized the xl residuals of the differences between the crossing times and lAE. There 
are two major error sources in the estimation of the timing differences. First, the calculated 
arm crossing times depend sensitively on the placement of the spiral arms, and we made a 
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comparison between the crossing times for our adopted model and that of Russeil (2003) to 
estimate the timing error related to uncertainties in the position of the spiral arms (approx¬ 
imately ±8 Myr except in the case of the crossing of the Scutum-Crux arm on the far side 
of the Galaxy where the difference is ~ 40 Myr). Secondly, there are errors associated with 
the estimated mid-times of the lAE, and we used the scatter between the various estimates 
in columns 2 - 5 of Table 1 to set this error (approximately ±14 Myr). We adopted the 
quadratic sum of these two errors in evaluating the xl statistic of each £t. The results of 
the htting procedure for various model and sample assumptions are listed in Table 2. 

The first trial £t was made by hnding the xl minimum that best matched the crossing 
times with the lAE midpoints from Shaviv (2003) (given in column 5 of Table 1 and noted as 
“Midpoint” in column 2 of Table 2). All four arm crossings were included in the calculation 
(indicated as 1 - 4 in column 3 of Table 2 ) that used the adopted model for the Galactic 
potential with a Galactocentric distance Ro = 8.0 kpc and and a stellar disk exponential scale 
length of Rd* = 2.4 kpc (model 7)^2 from Dehnen & Binney 1998b; see columns 4 and 5 of 
Table 2). The best £t difference (column 6 of Table 2) is obtained with — flp = 12.3 ±0.8 
km s“^ kpc“^, where the error was estimated by hnding the limits for which xl increased 
by 1. This £t gave reasonable agreement between the lAE and crossing times for all but 
the most recent crossing of the Sagittarius - Garina arm. Thus, we made a second £t {4^2 
in Table 2) using only the crossings associated with lAE 2-4, and this solution (with 
fl© — Qp = 11.9 ± 0.7 km s“^ kpc“^) is the one illustrated in Figure 3. The crossing times 
(given in the hnal column of Table 1) agree well with the adopted lAE midpoints. Our 
results are similar to the estimate of Oq — Qp= 10.4 ± 1.5 km s“^ kpc“^ from Shaviv (2003) 
who assumed a circular orbit for the Sun in the Galaxy. 

We also computed orbits using two other models for the Galactic potential from Dehnen 
& Binney (1998b) and determined the best £t spiral speeds for these as well. Fit ^3 in Table 2 
was made assuming a larger Galactocentric distance Ro = 8.5 kpc but with the same ratio 
of Rd-k/Ro (model 41^2h in Dehnen & Binney 1998b), and the resulting best £t spiral speed 
is the same within errors as that for our adopted model. We also computed an orbit for a 
potential with a larger value of disk exponential scale length Rdi,/Ro (model ^(^3 in Dehnen 
& Binney 1998b), but again the best fit spiral speed (£t in Table 2) is the same within 
errors as that for our adopted model. Thus, the details of the adopted Galactic potential 
model have little inhuence on the derived spiral pattern speed needed to match the lAE 
times. 

We might expect that the lAE midpoint occurs somewhat after the central crossing 
of the arm. For example, Shaviv (2003) suggests that the lAE midpoint may occur some 
21 - 35 Myr after the central arm crossing due to the difference in the stellar and pattern 
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speeds (so that the cosmic rays move ahead of arms as the stellar population does) and 
to the time delay between stellar birth and supernova explosion of the SN II cosmic ray 
sources. Furthermore, if ice ages are triggered by encounters with dense clouds as suggested 
by Yeghikyan & Fahr (2004), then the ice age may not begin until the Sun reaches the gas 
density maximum at the center of the arm. Thus, we calculated a second set of best fit spiral 
speeds to match the mean crossing and icehouse starting times (Frakes et ah 1992), and these 
are listed as fits ^(^5 and 7)^6 in Table 2. This assumption leads to somewhat smaller values 
of IIq — Qp, but ones that agree within errors with all the other estimates. 

We offer a few cautionary notes about possible systematic errors in this analysis. First, 
the £t of the lAE and arm crossing times depends on the difference hi© — Qp, and if our 
assumed value of IIq eventually needs revision, then so too will the spiral pattern speed hip 
need adjustment. For example, Reid & Brunthaler (2004) derive an angular rotation speed of 
^LSR = 29.5 ±1.9 km s“^ kpc“^ for the local standard of rest based upon Very Long Baseline 
Array observations of the proper motion of Sgr A* with respect to two extragalactic radio 
sources. If we suppose the local Galactic rotation curve is flat, then Gq = Qlsr 8.0/i?^ = 
28.7 ± 1.8 km s“^ kpc“^, where Rg = 8.23 kpc is the Sun’s mean Galactocentric distance. 
Adopting this value results in a spiral pattern speed of Gp = 16.8 ± 2.0. Second, our 
calculation ignores any orbital perturbations caused by close encounters with giant molecular 
clouds that cause an increase in the Sun’s motion with respect to a circularly rotating frame 
of reference. Nordstrom et ah (2004) present of a study of the ages and velocities of Galactic 
disk stars that indicates a net increase in the random component of motion proportional 
to time raised to the exponent 0.34. Thus, we would expect that the Sun’s random speed 
has increased through encounters by only Ri 4% over the last 500 Myr, too small to change 
the orbit or the arm crossing times estimates significantly. Third, we have ignored the 
deviations in the gravitational potential caused by the arms themselves. The Sun presumably 
slows somewhat during the arm crossings so that the duration of the passage is longer than 
indicated in our model, but since our model of the gravitational potential represents an 
azimuthal average, the derived orbital period and interarm crossing times should be reliable. 

Leitch & Vasisht (1998) argue that mass extinctions may also preferentially occur dur¬ 
ing spiral arm crossings. However, they proposed that a spiral pattern speed of Qp = 19 
km s“^ kpc“^ is required to find consistency between times of mass extinctions and spiral 
arm crossings, and if correct, then the relationship between ice ages and arm crossings would 
apparently be ruled out because Qp = 19 km s“^ kpc“^ is too large for the inter-arm crossing 
time to match the intervals between lAE (see Fig. 2 and Fig. 3). We show the times of 
the five major mass extinctions as X signs in Figures 1-3 (Raup & Sepkoski 1986; Benton 
1995; Matsumoto & Kubotani 1996). We see that in fact the lower value of Gp = 14.4 
km s“^ kpc“^ (Gq — Qp = 11.9 km s“^ kpc“^, as shown in Fig. 3) also leads to a distribution 



of mass extinction times that fall close to or within a spiral arm passage, so the association 
of mass extinctions with arm crossings may also be viable in models with pattern speeds 
that are consistent with the ice age predictions. 

Our calculation of the Sun’s motion in the Galaxy appears to be consistent with the 
suggestion that ice age epochs occur around the times of spiral arm passages as long as the 
spiral pattern speed is close to Op = 14 — 17 km s“^ kpc“^. However, this value is somewhat 
slower than the 20 ± 5 km s“^ kpc“^ preferred in recent dynamical models of the Galaxy 
(Amaral & Lepine 1997; Bissantz et ah 2003; Martos et al. 2004). The resolution of this 
dilemma may require more advanced dynamical models that can accommodate differences 
between pattern speeds in the inner and outer parts of the Galaxy (for example, a possible 
resonance between the four-armed spiral pattern moving with Op = 15 km s“^ kpc“^ and a 
“two-armed” inner bar moving with Op = 60 km s“^ kpc“^; Bissantz et al. 2003). 
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Fig. 1.— The Sun’s position in the Galaxy over the last 500 Myr expressed in cylindrical 
coordinates, R the distance from Galactic center {top), 0 the azimuthal position in the disk 
relative to 0 = 0° at present {middle), and Z the distance from the plane {bottom). Thick 
line portions mark icehouse epochs on Earth (Frakes et al. 1992), and X signs indicate times 
of large mass extinctions on Earth. The names of the geological eras and periods over this 
time span are noted at top. 
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Fig. 2.— A depiction of the spiral arm pattern of the Galaxy as viewed from above the 
plane. The plus sign marks the center of the Galaxy while the main four arms plus the local 
(Orion) spur are indicated as gray shaded regions. The dotted line through the center of the 
Galaxy indicates the location of the central bar (Bissantz et al. 2003). The Sun’s path in the 
reference frame of the spiral arms is indicated with a solid line (for Qp = 20 km s“^ kpc“^), 
and diamonds mark time intervals of 100 Myr back in time from the present {top diamond). 
The thick portions correspond to icehouse times and the X signs indicate times of large mass 
extinctions. 
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Fig. 3.— A depiction of the Sun’s motion relative to the spiral arm pattern in the same 
format as Fig. 2, but this time for a smaller spiral pattern speed (f2„ = 14.4 km s“^ kpc“^). 
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Table 1. Mid-points of Ice Age Epochs 


Ice Age Crowell (1999) Frakes et al. (1992) Veizer et al. (2000) Shaviv (2003) Arm Crossing (Fit #2) 
Epoch (Myr BP) (Myr BP) (Myr BP) (Myr BP) (Myr BP) 

1 . <22 < 28 30 20 80 

2 . 155 144 180 160 156 

3 . 319 293 310 310 310 

4 . 437 440 450 446 446 


Table 2. Fits of Spiral Arm Pattern Speed 


Fit lAE lAE Ro Rd* fl© — flp 

Number Times Sample (kpc) (kpc) (km s“^ kpc“^) 

1 . Midpoint 1 - 4 8.0 2.40 12.3 ± 0.8 

2 . Midpoint 2 - 4 8.0 2.40 11.9 ±0.7 

3 . Midpoint 2 - 4 8.5 2.55 11.8 ±0.6 

4 . Midpoint 2 - 4 8.0 2.80 11.8 ±0.7 

5 . Starting 1 - 4 8.0 2.40 11.6 ±0.8 

6 . Starting 2 - 4 8.0 2.40 11.4 ±0.6 





















